Glassy aspects of melting dynamics 

On melting dynamics and the glass transition, Part I 
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The following properties are in the present literature associated with the behaviour of super-cooled 
glass-forming liquids: faster than exponential growth of the relaxation time, dynamical hetero- 
geneities, growing point-to-set correlation length, crossover from mean field behaviour to activated 
dynamics. In this paper we argue that these properties are also present in a much simpler situation, 
namely the melting of the bulk of an ordered phase beyond a first order phase transition point. This 
is a promising path towards a better theoretical, numerical and experimental understanding of the 
above phenomena and of the physics of super-cooled liquids. We discuss in detail the analogies and 
the differences between the glass and the bulk melting transitions. 

PACS numbers: 64.70.Q-,05.50.+q,64.70.dj 



Almost any liquid becomes a glass when cooled fast 
enough ■ Many different scenarios and theories have 
been proposed over the time to describe the nature of 
glasses. Yet, it is still not known what is the fundamen- 
tal principle behind the experimentally observed abrupt 
change in the relaxation time of super-cooled liquids. Is 
there an underlying critical phenomenon behind the glass 
transition or not? Is the glass transition a thermody- 
namic or purely dynamic notion? What is the correct 
theory of super-cooled liquids? These questions remain 
unanswered and widely discussed. 

It is the experimental and numerical observations on 
super-cooled liquids that remind us of thermodynamic 
phase transitions. The most remarkable experimental 
facts about the transition from liquids to glasses are in- 
deed: (a) The extremely fast rise of the relaxation time 
T, that increases easily by several orders of magnitude 
as the temperature is decreased by a few percent. It 
is well approximated by the Vogel-Fulcher-Talman law 
T cx exp [A/{T - Tvft)] M- (b) The extrapolated tem- 
perature TvFT is found to be very close to the Kauzmann 
temperature Tk where the extrapolated entropy of the 
supercooled liquid becomes smaller than that of the crys- 
tal H. (c) The Adam-Gibbs relation r cx exp [C/AS'(r)] 
0, where AS'(T) is the difference between the two en- 
tropies, is observed with a good accuracy. 

It has been suggested that this quasi-universal be- 
haviour should be related to the existence of a length 
scale growing as the glass transition is approached, and 
the hunt for such a length has been a central theme in 
the field in the last decade. A purely dynamical quan- 
tity, proposed originally in the context of mean-field spin 
glasses [7j, uses a four-point density correlator in both 
time and space, and has led to the notion of dynamical 
heterogeneities and to the so-called dynamical suscepti- 
bility (usually refereed to as xa)- Indeed, the growing 
of dynamical heterogeneities has been observed in glass 
formers [1, Q. By contrast, other groups have consid- 



ered the possibility of a static (thermodynamic) growing 
length scale and proposed the point-to-set correlation, 
that is the correlation of a sub-system with its frozen 
boundaries [lol - [l^ . Recent numerical and experimental 
works confirmed that indeed there seems to be such a 
grow ing thermodynamic length scale in supercooled liq- 
uids [l3| . 

One of the standard routes to the above phenomenol- 
ogy goes through the Random First Order Theory 
(RFOT) (l3 - [T8l | according to which a good starting point 
of the glass phenomenology are mean-field s pin glasses, 
and in particular the p-spin glass model |19l422l |. Mean- 
field spin glasses have an interesting phenomenology: 
they behave like a liquid/paramagnet at large temper- 
ature, and when cooled down the equilibration time di- 
verges as a power law at a temperature Tmct and the 
relaxation process is described by the Mode-Coupling 
Theory (MCT) [23|-|28t. The true thermodynamic glass 
transition, however, arises at Tk < Jmct in these mod- 
els. According to RFOT, the correct theory of the glass 
transition is the finite dimensional counterpart of this be- 
haviour. A phenomenology called the the mosaic picture 
[6l[l0l. [T5l - IT7j is used in order to explain that in structural 
glasses the relaxation time diverges only at Tk- This 
crossover from the power-law divergence at Tmct (in the 
mean-field) to the super-exponential divergence at Tk (in 
finite dimension) is at the root of the RFOT theory. The 
validity of this picture is still, however, disputed (see for 
instance [l^, Ho] and references therein) . 

In this work, we argue that the above phenomena, usu- 
ally associated in the present literature with the dynam- 
ics of supercooled liquids close to the glass transition, 
also appear in a much simpler problem - the bulk melt- 
ing process of the fully ordered phase above an ordinary 
first order phase transition. All the following ingredi- 
ents arise: the crossover from a power-law divergence 
of the relaxation time to an activated dynamics and a 
Vogel-Fulcher-Talman-like divergence; the presence of a 



2 



growing length scale associated with point-to-set correla- 
tions; the divergence of the dynamical susceptibility and 
the presence of dynamical heterogeneities; the presence 
of a plateau in the dynamical correlation function with 
increasing life-time as the transition is approached. 

Analogies between the glass transition and the first or- 
der phase transitions are often evoked ~ as testified by 
the very name of the random first order theory. We ar- 
gue that the extend to which these analogies are valid 
and useful is larger than previously anticipated. The 
existence of this correspondance calls for more detailed 
theoretical, numerical and experimental investigations in 
systems with a first order phase transition. It offers a 
simpler way to understand some of the aspects and phe- 
nomenology of glassy dynamics, and inversly provides 
new ways to look at the bulk melting problem as well. 

In a subsequent paper (sTj , we shall go beyond a simple 
analogy and show that in some systems there is an exact 
mapping between the equilibrium glassy dynamics and 
the melting of ordered phase. This is true in particular 
in the mean fie ld jo-sp in models [l9l - l2]| . and also on the 
Nishimori line [3 21 133| in finite dimensional systems. 

This paper is organized as follows: In the first section, 
we briefly introduce the melting problem in systems with 
a first order phase transition, and stress the differences 
between surface melting and bulk melting. In the second 
section, we review the properties of bulk melting above 
a first order phase transitions on the mean-field level. 
As an example we use the exactly solvable p-spin fer- 
romagnet on a fully connected lattice, that undergoes a 
first order ferromagnetic phase transition. In the third 
section, we move to the finite dimensional case, and dis- 
cuss briefiy the nucleation and growth arguments to pre- 
dict properties of melting in finite-dimensional systems. 
In the fourth section, we simulate numerically the Potts 
model on a two dimensional lattice to show that indeed 
most of the phenomenology discussed in the glass transi- 
tion in finite dimensional systems arises generically also 
during the melting of the ordered phase. Section V ex- 
plains two crucial differences between dynamics of super- 
cooled liquids close to the glass transition and melting in 
a generic system with a first order phase transition. In 
the subsequent paper [31] we then study systems with a 
first order transition where these differences disappear. 
We finally review our results and discuss some criticisms 
of the RFOT scenario in the light of our findings. 



I. BULK MELTING 

In this section we specify what we mean by melting 
process above a first order phase transition, and in par- 
ticular we emphasize that we deal here with bulk melting 
[U, instead of the much studied — and arguably more 
practically important — surface melting process. 

Consider a system with a first order phase transition, 
for instance the solid-liquid transition in crystalline struc- 
tures or in three-dimensional hard spheres [SS^, or any 



given ferromagnetic spin system with a first order tran- 
sition such as the Ising model in field or the Potts model 
in temperature [s^ . Start in the fully ordered state (the 
crystal for solids, the perfect packing for hard spheres, 
or the completely magnetized state in the spin models) 
and suddenly change the pressure/temperature/field in 
order to put the system in the liquid/paramagnetic/less 
ordered phase. The system will melt. How long does this 
melting process take, and what are its properties? How 
does a solid turn into a liquid! 

The occurrence of a discontinuous transition suggests 
that melting of crystalline or solid matter in nature is 
a first order phase change that requires a finite latent 
heat, and a nucleation mechanism. This is indeed the 
case provided some precautions are taken. Consider for 
instances solids, unlike super-cooling of liquids, super- 
heating of crystalline solids is observed to be extremely 
difficult if not impossible [33.l37l.[3^ in most experiments. 
This peculiarity is due to the fact that melting of a crys- 
tal kept at a homogeneous temperature always begins on 
its free surface. Surface melting is indeed the dominant 
mechanism for melting of solids in nature, as it is a much 
faster process than the bulk melting. It was, however, 
realized that by suppressing surface melting [4l| super- 
heating to temperatures well above the equilibrium melt- 
ing point could be achieved. In this case, one recovers 
the nucleation phenomenology of first order phase tran- 
sitions, and this is the kind of melting processes that, 
as we will show, displays many analogies with glassy dy- 
namics. Such bulk melting problems and the study of the 
super-heated solids have received a boost of experimental 
and theoretical studies recently [4l|, . 

Here, we shall work with very simple models, far from 
realistic atomic or sphere systems. Our motivation stems 
from the fact that we expect properties of bulk melting to 
be universal and to be qualitatively similar in solids, mag- 
netic systems, hard spheres or polymers with a first or- 
der transition. Following the classical mapping between 
spin systems and lattice gas [43|, we will concentrate on 
the first order transition in Ising and Potts spin models. 
Since these models are more easily amenable to analytic 
studies and simulations, this will allow us to discuss uni- 
versal behaviour, such as spinodal and nucleation mech- 
anisms. Moreover, in theoretical considerations surface 
melting can be suppressed by simply using fixed or pe- 
riodic boundary conditions. Finally, many activities in 
studies of the glass transition have been devoted to spin 
models, this is hence a natural setting for discussing the 
analogies between glassy and melting dynamics. 



II. MELTING IN MEAN-FIELD SYSTEMS 

In this section, we start by discussing the bulk melting 
transition at the mean-field level using a simple solvable 
model. An elementary example of a mean-field system 
with a first order phase transition is the fully-connected 
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FIG. 1: (color online): Gibbs free energy /(m, /?) as a func- 
tion of the magnetization m for different temperatures /3 for 
the ferromagnetic fully-connected 3-spin model. The figure in- 
cludes the free energy at the critical temperature /3c = 0.67209 
(blue), at the spinodal temperature (3s = 0.57217 (red), and 
at the temperature where the free energy starts to be a convex 
function of the magnetization /?* = 0.435 (green). This illus- 
trates the generic behaviour of the free energy in a mean-field 
system with a first order phase transition. 



ferromagnetic p-spin Ising model, for p = 3 

n^-JJ2s^s,Sk, (1) 

ijk 

where the sum is over all values of indices i,j,k. This 
is a simple generalization of the Curie- Weiss model. 



A. The static behaviour 
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FIG. 2: (color online) The melting process starting from the 
completely magnetized m(0) = 1 state for different temper- 
atures in the fully connected 3-spin ferromagnet. Bottom: 
The time evolution of the magnetization m{t), from left to 
right, 13 = 0.35, 0.4, 0.45, 0.5, 0.52, 0.54, 0.56, 0.565, 0.57, 
0.571, 0.572, 0.5721 and the spinodal temperature (in red) 
Ps = 0.57217. m{t) shows a two-steps relaxation when ap- 
proaching the spinodal point (that starts to be observable 
at roughly /3«). The time needed for the system to melt di- 
verges according to eq. ((SJ (see inset). Top: The dynamical 
ferromagnetic susceptibility eq. ([7} for the same set of tem- 
peratures. It has a maximum around r and converges at large 
times to the static value xf = P- The peak of the suscepti- 
bility diverges when approaching the spinodal point as eq. ^ 
(see inset), indicating the growth of a correlated volume. 



Just like the Curie- Weiss model, the ferromagnetic p- 
spin model is exactly solvable on a fully connected lattice 
(which is a very crude mean-field approximation to the 
finite dimensional case). To ensure extensivity of the 
thermodynamic potentials we take J = 1/NP~^. The 
Gibbs free energy as a function of magnetization m and 
inverse temperature /? reads (see the appendix): 

/3/(m, (3) — —fimF — log [2 cosh atanh(m)] + m atanh(m) 

(2) 

while the self-consistent equation for magnetization reads 

m = tanh (/3pm^~^) . (3) 

As illustrated in Fig. [Tl for p = 3, this yields a first 
order phase transition at /3c = 0.67209. The low temper- 
ature ferromagnetic phase is, however, locally stable up 
to the spinodal point at /3s — 0.57217. The Gibbs free 
energy f{m,l3) starts to be a convex function of magne- 
tization above temperature corresponding to /3* = 0.435. 



B. The dynamical behaviour 

Consider the following Glauber dynamics: At each in- 
terval of time we pick up one spin, and set it -1-1 
with probability p+ = (l-|-tanh/3/i)/2, and —1 with prob- 
ability p- = {1 — tanh/3/i)/2, where h = p-wP~^ is the 
local field on the spin. This dynamics is also exactly solv- 
able for the fully connected ferromagnetic p-spin model 
(see again the appendix for derivations). The average 
(over realizations) of the magnetization {m{t)) evolves in 
time according to the following differential equation 

^^^^^^ = -{m{t))+ tanh [pp{m{t))P-^] . (4) 

Fig. [2] shows how the magnetization evolves in time if the 
system is initialized in the ferromagnetic configuration 
m{t = 0) = 1 for different temperatures. As the spinodal 
temperature of the ferromagnetic phase Ps is approached 
from above, we observe a two-steps relaxation. The mag- 
netization of the flat region is roughly the magnetization 
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of the appearing metastable state and its length defines 
the relaxation time. This relaxation time diverges as a 
power law 



r oc 



(5) 



The flattening of the relaxation process starts roughly at 
the temperature /3* . The growing plateau in the decay of 
the order parameter and the power-law divergence of the 
relaxation time are classical features of a melting process 
above a mean- field first order transition H, |43 . 



C. Static and dynamic diverging length scales 

Coming back to the discussion in the introduction, one 
may ask if there is a growing length, or volume, associ- 
ated with the divergence of this relaxation time. The 
naive answer is no: The magnetic susceptibility that di- 
verges at a second order phase transition stays indeed 
finite in systems with a first order phase transition when 
the ferromagnetic spinodal temperature Ts or the criti- 
cal temperature Tc is approached from above. This an- 
swer is, however, wrong. Consider the following (time- 
dependent) correlation function, averaged over many re- 
alizations of the melting dynamics 



a,{t) ^ {s.{t)s,{t)) - {s.{t)){s,{t)) 



(6) 



It is natural to ask how many spins have correlated moves 
at a given time t. An estimates of this is given by inte- 
grating this function over the whole system 

One recognizes that, for large time, this is nothing but 
the equilibrium magnetic susceptibility (we have used the 
multiplication by /3 in order to be consistent with the 
usual definition of the magnetic susceptibility), and we 
shall thus call this quantity the dynamic magnetic sus- 
ceptibility. Certainly the equilibrium magnetic suscepti- 
bility is finite in a system with a first order transition 
when the ferromagnetic spinodal temperature is ap- 
proached from above. Let us look, however, to the finite 
time behaviour of XF{t), it follows a differential equa- 
tion derived in the appendix. The solution for different 
temperatures is plotted in Fig. [5] 

We see that XF{t,(3) increases with time, reaches a 
maximum roughly at the relaxation time t = T{f3), and 
then decays to reach its equilibrium value (which is noth- 
ing else than the Curie law for a paramagnet xf = /?)■ 
The maximum value of xf{t,I3) grows, however, very 
fast as the spinodal point /3s is approached. The num- 
ber of spins that are correlated in the melting process is 
in this system maximum roughly at the relaxation time 
t — t{I3), and thus diverges at the spinodal point as 



Xf{t,(3) cx ( 1 - ^ 



This defines the dynamic length (or volume) that diverges 
at the spinodal in the melting problem. The origin of the 
critical exponent 5/2 is elucidated in a recent work [46j . 

Is there also static volume/length diverging at /3s? 
Again, the answer is positive. In order to see it, we con- 
sider the system with boundary conditions fixed to pos- 
itive magnetization. In the paramagnetic phase above 
the spinodal temperature, fixing the boundaries in this 
way does not inffuence the equilibrium bulk properties. 
But how far should the boundaries be to find back the 
zero magnetization phase? This is a well defined point- 
to-set correlation length. In a fully connected model it 
is, however, rather tricky to define boundary conditions. 
Let us thus consider the model on a tree of fixed finite 
degree. The model can then be solved with the Bethe- 
Peierls method (see the appendix). In the large degree 
limit the model on tree becomes equivalent to the fully 
connected one, and the magnetization at distance I from 
the boundaries follows the recursion 



mi+i 



tanh 



(9) 



In Fig. [3] we show how the magnetization evolves with 
distance from the boundaries. One observes that the dis- 
tances one has to go in order to recover the zero magne- 
tization diverges when approaching the spinodal point. 
Actually, since eq. © is simply the discrete version of 
eq. ([1]), this length diverges as 



1 - 



(10) 



(8) 



The length £ we have just defined is in fact well know 
in the theory of wetting problems (47l |. where indeed a 
divergence is expected. Let us point out that the very 
same diverging length on a tree can be defined in an- 
other equivalent way: One takes a finite tree of depth 
^ and computes how large ^ needs to be such that the 
magnetization in the center (or the total magnetization) 
is bellow a given threshold. On a tree the two definitions 
are equivalent and the length £ is equal to the length 
^ as the two processes are described by the very same 
equation. However, the two definitions differ in finite di- 
mensional systems and we will see that it is the second 
length ^ that is associated to the point-to-set correlation 
length in finite dimension (while the first one £ is the 
finite dimensional wetting length). 

Using the fully connected ferromagnet, we have here 
described the basic features of the melting process in a 
mean-field setting. The time r needed to decorrelate di- 
verges as the temperature approaches the spinodal point. 
This translates into a growing plateau in the decay of 
the order parameter. The divergence of a time scale is 
accompanied by at least two diverging length scales: the 
dynamical one, that shows how many spins are corre- 
lated in the melting process at the relaxation time scale 
T, and a static one, that corresponds to the way the sys- 
tem decorrelates with ordered boundary conditions. 
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FIG. 3: (color online) The magnetization m; at dis- 
tance I from boundaries fixed to positive magneti- 
zation in the paramagnetic phase of the mean-field 
ferromagnet. From left to right, we show j3 — 
0.5, 0.51, 0.52, 0.53, 0.54, 0.55, 0.56, 0.57, 0.571, 0.572, 0.5721 
and (in red) /3s — 0.57217. The distance one has to go from 
the boundaries in order to recover the equilibrium magne- 
tization diverges as eq. (fTU)) as the spinodal is approached. 



The learned reader will recognize that these two 
lengths are nothing but the dynamical four-point suscep- 
tibility and the point-to-set correlation in disguise. This 
is the message we want to convey: these length scales de- 
fined initially in the mean field theory of the glass transi- 
tion do also diverge in the simpler situation of the melting 
of the ordered phase above a first order transition. 

At this point, a word of caution should be made: In 
order to observe these diverging length and time scales 
when lowering the temperature down to the ferromag- 
netic spinodal Ts, it is necessary to start from the or- 
dered phase and to consider its melting towards the high 
temperature phase. If one instead starts from the equilib- 
rium liquid/paramagnet state, then the phenomenology 
is much less interesting: the equilibrium relaxation time 
in the paramagnetic phase does not show any particular 
interesting feature, apart from increasing when approach- 
ing the eventual low temperature paramagnetic spinodal 
point beyond which the system orders (crystallizes). Of 
course, if this ordering is avoided a true glass transition 
could be observed at yet lower temperatures, but this is 
a different phenomena from the one described above. 



III. BASICS OF NUCLEATION THEORY 

In this section, we discuss what happens in finite di- 
mensional systems and recall briefly the classical nucle- 
ation theory in order to explain that all the diverging 
time and length scales discussed above should diverge at 
the critical point in finite dimensional systems (and not 
anymore at the spinodal - which is no longer thermo- 



dynamically well defined). In the next section, we will 
illustrate these properties for the two-dimensional Potts 
ferromagnet. 

In finite dimensional systems the free energy density is 
always a convex function of the magnetization, hence the 
life-time of any metastable state is always finite. Phys- 
ically, this is due to the nucleation of the equilibrium 
phase that is absent in mean-field geometries. Metasta- 
bility hence becomes a purely dynamical notion and for 
this very reason giving a proper theoretical definition of a 
metastable state in finite dimension can be tricky 48, 49| . 

Nucleation theory, introduced originally by Gibbs [50j , 
describes the way metastable states decay beyond a first 
order critical point. The metastable phase has higher free 
energy, hence it must decay by nucleation of droplets of 
the stable (lower free energy) phase. Nucleation results 
from a competition between the bulk free energy differ- 
ence and the surface tension between the two phases. 
When a system is in a metastable phase the free-energy 
cost AF{£) of a droplet of size £ is given by 



AF{£) = -vj^sf) + sd'^-^r ., 



(11) 



where Vd and Sd are the volume and surface factors re- 
lated to the shape of the droplet in dimension d, Sf is the 
free-energy difference between the two different phases, 
and r is the surface tension. The free energy cost is 
maximum at size 



/ — 



id-i)Sd r 

dVd Sf 



(12) 



Droplets smaller than ^max have the tendency to shrink, 
while those larger than ^max will expand fast (the later 
process is often described by the Kolmogorov-Johnson- 
Mehl-Avrami theory [5l|). In order to relax to equi- 
librium, thermal fluctuations must cross the free ener- 
getic barrier A_F'(^max)- This nucleation time is often 
estimated by the Arrhenius formula 



« exp[/3AF(£ 



exp 



(13) 



where Ad is a dimension dependent constant. The first 
order melting point is defined by (5/ = and therefore 
both the nucleation time (jl3l) and the nucleus size 
diverge at the first-order critical point. 

When a critical droplet is formed, it will grow and in- 
vade the system. The relaxation process thus depends 
on both the nucleation time and on the growth of the 
nucleus. Moreover relaxation can follow from the ap- 
pearance of many droplets in the system: Indeed the 
nucleation time t^ucI gives the average time needed to 
nucleate a critical droplet at a given place in the system. 
If, however, the system is much larger than the critical 
nucleus size (a condition that becomes harder and harder 
to fulfill as we approach the first-order transition point), 
the relaxation process will not be done via the growth 
of a single droplet but instead will result from the nucle- 
ation of many droplets that will simultaneously grow and 
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invade the system. Whether the first or the second sce- 
nario is happening depends on whether the system size is 
much larger than the critical nucleus size. In both case, 
however, the growth time has to be added to the nucle- 
ation time to obtain the total relaxation time needed to 
exit from the metastable state. 

Consider a large system, much larger than the size of 
a critical nuclei. Given the probability to nucleate a crit- 
ical droplet by a unit of volume and time is proportional 
to p cx e"'^^^^'"""') oc 1/rnuci, in a very large system 
there should therefore be roughly one droplet per vol- 
ume Tnuci after a finite time. In most (non disordered) 
systems, and for most dynamics, the droplet radius is 
expected to grow as a power law i{t) oc t^/^ with time. 

Therefore, after a time of order r^^^^j droplets will start 
to touch each other and the portion of the system in the 
new phase will percolate. This predicts that the relax- 



ation time is T^eiax oc exp 



i.e. it behaves just 



_d (SfY- 

as the nucleation time, up to a constant rescaling of the 
barrier. This is the first important conclusion: the relax- 
ation time diverges super-exponentially at the first order 
transition point. Note also that at this time the system 
is maximally heterogeneous as a finite portion is the new 
phase, while the rest is in the initial one, and we shall 
see that indeed dynamical heterogeneities are maximal 
at this point. 

This presentation is oversimplified, as many aspects of 
the nucleation and growth process are neglected (for in- 
stance the possibility to have non-compact droplets), and 
from a generic point of view nucleation theories could 
predict more complex exponents in eqs. (1121113^ . see for 
instance [44]. The bottom- line is, however, that when 
we consider the melting process in finite dimension, the 
power-law divergence at the spinodal point transforms 
into super-exponentially slow relaxation at the first order 
transition temperature Tc. Consequently, even before the 
transition the relaxation time exceeds the experimental 
time and metastability becomes observable and unavoid- 
able in practice. This behaviour of the relaxation time 
at a first order phase transition is thus extremely close to 
the Vogel-Fulcher phenomenology for glasses [j]. More- 
over, if we consider a first-order phase transition driven 
by entropy (as e.g. the 3D hard sphere problem) the free 
energy difference in eq. (jl3p is replaced by the entropy 
difference: in this case the formula looks just like the 
Adams-Gibbs relation for glasses 



IV. MELTING DYNAMICS IN THE 2D POTTS 
FERROMAGNET 

We shall now illustrate the above behaviour on the 2D 
ferromagnetic Potts model, which is one of the simplest 
models with a first order phase transition. It is defined 
by the Hamiltonian 



n = -Y,5{s^,,s,), 



(14) 



where the sum is over all neighboring spins in the two- 
dimensional square lattice, and where S = 1, 2, . . . , q. For 
q = 2, it reduces to the Ising model. In two dimensions, 
one can show [36*1 that the critical temperature separating 
the ordered phase from the paramagnetic one is 



Tc{q) 



1 



log (1 + Vg) ■ 



(15) 



The transition is a first-order one for q > A. In order not 
be too close to a second order phase transition, we shall 
work with q — 10. In this case, the critical temperature 
is Tc(lO) « 0.70123. At zero temperature, the ordered 
configurations are the 10 ones where all spins take the 
same value. We shall select one of them, e.g. all S'i = 1, 
and study the melting of this configuration at different 
temperatures. We simulate a lattice of size L — 500, 
i.e. N — — 250000 spins and use periodic boundary 
conditions (except for the point-to-set correlation). 



A. Melting and relaxation time 

We start by repeating the procedure we used in the 
mean-field ferromagnetic model: We initialize the system 
into a fully ordered configuration Si = 1 and study the 
dynamics at temperature T > Tc (for the simulations, 
we used the heat bath dynamics) . In Fig. |4] we show the 
time dependence of the normalized magnetization 



m{t) 



t),i 



N{q-1) 



(16) 



It develops a plateau when approaching the phase tran- 
sition, as typical in glassy systems. The plateau grows 
super-exponentially fast when T — >■ Tc, see the inset of 
Fig. m as predicted by the nucleation theory. 



B. Point-to-set like correlations 

We shall now fix the boundaries to be in the fully or- 
dered configuration and study the magnetization they 
induce for temperature T >Tc. We follow [l^ and con- 
sider a box of size L with fixed fully magnetized boundary 
conditions. We then run a Monte-Carlo simulation, us- 
ing the Wolff's cluster algorithm 52| in order to reach 
equilibrium [75J, and measure the total magnetization in- 
side the systems for different sizes L and temperatures T. 
The results are shown in Fig. [5j We see clearly a growing 
point-to-set correlation length that diverges as the transi- 
tion at Tc is approached. Defining the correlation length 
as the moment when the magnetization falls bellow the 
plateau (in practice, we used m < 0.8), we observed that 
the length grows approximately as l/{T—Tc)^-^^-^\ which 
is not very far from the linear scaling of the size of a 
critical nucleus Indeed, ^(T) is nothing but an es- 

timation of the nucleus average size [ssl ] in the melting 
process. As the temperature drops, melting dynamics is 
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t 



FIG. 4: (color online): Dynamical susceptibility (top) and 
decay of the magnetization m{t) (bottom) in the melting of 
the 2D 10-state Potts model with N = 250000 spins. From 
left to right: T = 0.8, 0.79, 0.78, 0.77, 0.76, 0.75, 0.74, 0.73, 
0.72, 0.71, 0.706, 0.704, 0.703, 0.7025 and 0.7022. The insets 
show the super-exponential divergence of the relajcation time 
(bottom) and of the maximum susceptibility (top), that is 

found to grow as ^max 0^ 



initiated by the reversal of larger and larger clusters of 
spins (the nuclei that will later expand), and the size of 
these clusters diverges at Tc- 

A word of caution is needed here: we have identified 
a diverging length scale in first-order transition process. 
But this length not be confused with the equilibrium cor- 
relation length, which does not diverge at a first order 
transition (in fact it reaches a value of w 7.5 lattice spac- 
ings exactly at the transition in this model, see [53] ) . The 
point-to-set length diverges while the equilibrium corre- 
lation one does not. 

To echo the discussion in sec. Ill C[ let us further em- 
phasize that the point-to-set construction we have used 
is also different from taking first a very large box with or- 
dered boundary conditions and then measuring the mag- 
netization at distance L from the boundaries. Indeed, 
for box sizes L ^ ^ nucleation is always favorable and 
the nucleus thus invades almost the whole system. From 
the data presented in Fig. [S]we can see that the magne- 
tization for L ^ ^ starts to decay as the inverse of the 
size L. This can be understood easily: If the correlation 
with the boundaries in a very large box is set by a char- 
acteristic distance £, then the local magnetization will be 
roughly one at distance d < i from the boundaries, and 
roughly zero at distance d > £. The total magnetiza- 
tion should thus behave as m cx £/L when L ^ ^. The 




20 40 60 80 100 



L 

FIG. 5: (color online) Equilibrium average magnetization in 
a box of size N — with fully magnetized fixed boundary 
conditions in the 2D Potts model. These data indicate the 
growth of the point-to-set length scale as the transition at 
Tc is approached. Defining the correlation length ^(T) such 
that m(i) < 0.8, we find that ^ oc (T - Tc)""'*'^', which is 
reasonably close to the prediction 1/(T — Tc) of nucleation 
theory . 



characteristic length £ is, however, much smaller and un- 
related to the point-to-set construction (in fact, it is in 
the so-called wetting length, which is expected, at least 
in simple mean field models, to diverge logarithmically 
at the transition, see for instance This is a major 

difference with the mean field case where at the spinodal 
point, both the correlation length with frozen boundaries 
in an infinite system (that is, the wetting length) and the 
length defined by the onset of magnetization in a finite 
system (that is, the point-to-set length) are equal and 
diverge in the very same way. This is not true anymore 
in finite dimensional systems. 



C. Dynamical heterogeneities 

We also computed the dynamical susceptibility XF{t), 
see Fig. S) Again, it yields the qualitative behaviour of 
the X4 observed in glassy systems, meaning that the dy- 
namical heterogeneities are the largest when the suscep- 
tibility reaches its maximum and they grow as the tem- 
perature is approaching the phase transition. As we are 
working with a 2D system, visualization is convenient 
and we thus depict these heterogeneities in Fig. [51 

Since the melting proceeds via nucleation, one might 
think that the maximum of the Xpit) denotes the vol- 
ume of the critical droplet ((T2|) . This is not the case, 
in fact the volume described by the maximum of xf is 
much larger than the critical nucleus size. This is clearly 
visible in the simulations depicted in Fig. [51 where the 
size of a typical critical nucleating droplet is about 10^, 
while the typical size of the correlated regions close to 
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FIG. 6: (color online) Nucleation process in the melting of 
the 2D Potts model with g = 10 at temperature T = 0.7025 
with A'^ — 1500'^ spins. The red dots are the spins that have 
the same value as the initial state S = 1, white are all the 
other spins. The times are from top to bottom 25000, 75000, 
120000, 220000, 260000, 380000. On the second image, we 
can already recognize some nucleating critical droplets that 
start to grow fast for larger times. The nucleating droplets 
percolate between images n. 4 and n. 5, this corresponds to 
the time at which the heterogeneities and the dynamical sus- 
ceptibility are the largest. The typical size of heterogeneities 
is much larger than the size of critical nuclei . 



the maximum of the xf is about 100^. 

The reason beyond this is simple and illustrates well 
the nucleation and growth mechanism. As described in 
Sec. HID after a time of ©(rioiax), grown droplets will 
touch each other and percolate; this is the moment at 
which the maximum of the xf arises, and we expect that 
this maximum thus behaves algebraically with the relax- 
ation (and the nucleation) time so that xf oc r", which 
is indeed what we find in the inset of Fig. |4l This shows 
that the nature of the static and the dynamic length scale 
is completely different, and that the dynamical suscepti- 
bility can probe much larger scales. The first one probes 
the length associated with critical nuclei while the second 
one depends on the growth process. Note that the expo- 



nential divergence of xf (with respect to the temperature 
difference as the critical point is approached) compares 
well with what is observed in many models with kineti- 
cally constraint or facilitated dynamics (55l - [57| . 

Given the simplicity of its origin, a law xf oc r" should 
indeed be quite generic for systems with nucleation and 
growth. However, the value of the exponent a depends 
on the type of growth and on the precise properties of the 
system. The growth process of nuclei can be in principle 
much slower in system with disorder (where, due to the 
pinning of the interface, activation is again necessary to 
grow the nuclei) with respect to the relatively fast process 
observed here, and its interplay with nucleation can be 
more complicated. In fact, the growth can probably be 
so slow that the two length scales might be comparable, 
and this will be investigated in a companion paper (3lj |. 



V. DIFFERENCES BETWEEN BULK MELTING 
AND GLASSY DYNAMICS 

So far we concentrated on the analogies between glassy 
dynamics and bulk melting, it is time to list several key 
differences. First of all, the analogy does not extend to 
non-equilibrium properties such as aging in glasses pGj . 
We are comparing here the melting process with the 
equilibrium-initialized dynamics in the glassy systems, 
but the out-of-equilibrium behaviour of glasses is very 
different. But let us concentrate on the differences in the 
dynamics from equilibrium. 

A crucial point is that there is no latent heat associ- 
ated to the glass transition, as the energy is continuous, 
while there is a latent heat in a general first-order transi- 
tion (for instance both in the ferromagnetic p-spin model 
and the 2D Potts model we discussed in this paper). To 
smear away this difference one should consider a first or- 
der transition driven purely by entropy, if the energy is 
continuous at the first order phase transition then there 
is no latent heat. This is the case for instance in a system 
of hard spheres. 

Another major difference in the dynamical behaviour 
is that the melting process happens once for all. When 
the system has melted, it stays in the liquid phase. How- 
ever, the equilibrium glassy dynamics is a stationary pro- 
cess which is time-translationally invariant. This is why 
the equivalent of the nucleation theory in glasses, the so- 
called mosaic picture [13, fl5l - [T7l | , is more complicated 
(and not yet fully understood). 

Our results seem, however, to point towards the glass 
transition being a melting of some special sort. In the 
glass problem, following the ideas of Goldstein [6l|, we 
can assume that we have a more complicated landscape, 
with many local minima. The glassy dynamics consists 
of a continuous and stationary melting of one minima to 
another one. This view is precisely the one adopted in 
the mosaic picture and in the RFOT [lol. IH - OtII . 

The comparison between melting and glassy dynamics 
will be pushed beyond a simple analogy in a compan- 
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ion paper [3l|, where we show that in a certain class 
of spin systems (both mean-field and finite dimensional 
ones) with a first order phase transition, one observes a 
purely entropic transition and the melting dynamics is 
a in fact equivalent to the equilibrium dynamics. In or- 
der to do so we explore properties on the Nishimori line 
[32,|3l,|6l,[6ll and study consequences of such a mapping 
for the physics of the glass transition. 



VI. DISCUSSION 

In this paper we explored analogies between the bulk 
melting dynamics above a first order phase transition 
and the equilibrium glassy dynamics. We showed that 
many properties associated with the dynamics of struc- 
tural glass formers are already present in the melting 
case. We have investigated several features of modern 
theory of super-cooled liquids — such as the diverging 
relaxation time, the plateau in the correlation function, 
dynamical heterogeneities related to the dynamical sus- 
ceptibility X4, the point-to-set correlations — and showed 
that they all can be easily recovered and understood 
in a much simpler setting of systems with a first order 
phase transition. In particular the dynamical — or mode- 
coupling — transition corresponds to the spinodal point, 
while the Kauzmann transition corresponds to the first 
order ferromagnetic phase transition. 

Noticing similarities between the glass transition and 
the first order one is of course not new. In fact, it is 
somehow implicit in the random first order theory of 
glasses [l5l - [l7| , where according to the mosaic picture ev- 
ery glassy state melts into another one and so on. The no- 
tion of nucleation stands at the roots of the theory deriv- 
ing the Adam-Gibbs law and the mosaic theory is nothing 
but the counter-part of nucleation theory for glassy tran- 
sitions. The construction of Franz and Parisi [59| and 
the earlier work of fTJ| are attempts to root the theory of 
the glass y tr ansition into a first-order setting. The recent 
work of [60| also accentuates the analogy, and the very 
definition of the point-to-set correlation was recently pro- 
posed as a tool to investigate nucleation in systems with 
a generic first order transition [53l | , and this is by far not 
an exhaustive list. The bottom-line of our approach here 
is that it is rather interesting to take this analogy liter- 
ally, to concentrate on the simpler melting problem, and 
to understand in detail its glassy aspects. We believe this 
approach can be interesting for the bulk melting problem 
per se. 

Before concluding, let us make a last comment. The 
random first order theory is not fully accepted as a good 
theory of the glass transition. Without taking a position 
on the question whether RFOT is the good theory for the 
glass transition or not, we want to point out here that 
many of the criticisms of RFOT apply equally to systems 
with an ordinary first order phase transition. Imagine for 
a moment that we would study melting dynamics with- 
out knowing about the underlying first order phase tran- 



sition. In this case researchers could be trying different 
fits to extrapolate the divergence of the melting time at 
a finite or zero temperature [g^. They could also criti- 
cize the use of ill-defined metastable states [l^, [s^ [6^ . 
Mean-field theories would predict a spinodal line with a 
power-law divergence (23l . [23 | , and the crossover to finite 
dimensional behaviour could be thought of as a myste- 
rious one (soj . There would be some claims that mean- 
field theories do not incorporate geometric properties and 
hence are not very relevant to finite dimensional systems 
[29I I . Some would concentrate on the heterogeneous fluc- 
tuations in the decay of the order parameter, and on the 
peculiar shape of critical nucleating droplets as the cen- 
tral object of interest [6^; claiming at the same time that 
there is no need to think about thermodynamics \6l\ and 
that melting should be modeled with purely dynamic sys- 
tems with simple rules [68j . In cases of systems such as 
diamond the melting transition might even be not con- 
sidered at all because the true equilibrium state is simple 
carbon — just as the crystal is always the true equilib- 
rium state in glass forming systems Luckily, first 
order phase transitions are much more established (since 
in this case, contrary to what happens in glasses, we can 
equilibrate at low temperatures), hence the answers to 
the above arguments are mostly well understood. Still, 
this shows that starting from a mean-field analysis and 
correcting it with nucleation-like arguments, as is done 
in RFOT, is not such a bad strategy a priori. 

Finally, a lesson from this study is that one should 
maybe look back to the first order phase transitions with 
the eyes of glassy phenomenology, both from an experi- 
mental and theoretical point of view. The following ques- 
tions are for instance appealing: How do the static and 
dynamic lengths behave in other models with first or- 
der transitions (with or without disorder)? What does 
the mode-coupling approximation predict for a first or- 
der phase transition? Can all the recent experimental in- 
vestigations of static and dynamic length scales in glass 
formers be repeated in a first order melting process? The 
way a solid fiows under shear is clearly related to the 
melting mechanism (H^, does the dynamical susceptibil- 
ity diverge as a power of the relaxation time as well? 

These questions are particularly appealiiig in the con- 
text of bulk melting in super-heated solids [41j . It would 
be for instance interesting to apply the dynamical corre- 
lation and susceptibility analysis of sections fll Cl andl lV CI 
in order to identify dynamical heterogeneities (which 
seem indeed to be present in this problem, see for in- 
stance 42] ) . This could also help to understand the inter- 
play between nucleation and growth in melting processes. 
Other common behaviour can be identified upon inspec- 
tion; compare for instance the string-like motion [70l | 
observed in glassy dynamics with the highly correlated 
ring and loop atomic motions behaviour observed in bulk 
melting (7]| (and more generally with the stringy nuclei 
observed close to pseudo-spinodal during first-order tran- 
sition [Z^). Another example is given by the behaviour 
of the shear modulus: according to a recent work within 
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the RFOT, it is predicted to vanish continuously at the 
mode-couphng transition point in glasses [7^, a state- 
ment clearly reminiscent of the Born criterion [t^I that 
defines the spinodal point of superheated solid as the 
moment where shear modulus is going to zero. Finally, 
studying the system with respect to frozen boundary con- 
ditions as in section IIVBI should also be interesting in 
order to estimate nucleation barriers numerically in the 
spirit of [5^, and to investigate the much studied prob- 
lem of the limit of the super-heated state. 

To conclude, there is a lot to learn on the crossroad be- 
tween glasses and first order melting, and we hope our ar- 
ticle will trigger new works in this direction. This will be 
further accentuated in a companion article where we pro- 
vide a set of models where the melting process is exactly 
equivalent to its equilibrium dynamics |3l| and where we 
show that the mean-field theory of the glass transition is 
mappable to a (special) melting problem. 



The self-consistent equation on m thus reads 

m = tanh(/3pTO^'~^ +/3/1) . (19) 

We are interested in the free energy as a function of the 
equilibrium magnetization m. This is obtained by the 
Legendre transform of eq. pS]) . that is, by /(m,/3) = 
5(/i*, /3, to) -|- h*m, where h* is given by the condition 
(fT9l) . It yields finally 

/3f{m,(3) = — /3m^— log [2 cosh (atanhTO)]-fTO atanh(TO) . 

(20) 

This formula was used to produce Fig. [TJ For p = 3, 
this solution yields first-order phase transition at /3c = 
0.67209. The low temperature phase, however, is locally 
stable until the spinodal point at /3s = 0.57217. 

2. The Bethe-Peierls solution 
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Appendix: The ferromagnetic p-spin model 

In this appendix, we remind how to solve the static 
and the dynamic behaviour of the fully connected ferro- 
magnetic p-spin model defined by Hamiltonian ([1]) with 
J = 1/7VP-1. 



1. The thermodynamic solution 



Consider the p-spin model on a tree with coordination 
number c. In this case, one can write an iterative re- 
cursion using the Bethe-Peierls strategy that relates the 
magnetization at level I + 1 with the one at level I [2l| 



tanh I (c — l)arctanh tanh (/? J)to^ ^ |. (21) 



Let us now take the limit of a very large connectivity and 
use J = p/c. In this limit 



mi+i 



tanh 



^p/?TO^ ^ 



(22) 



The stationarity condition of this equation leads to 
eq. (|19p. This formulation allows to define a correlation 
length as is done in sec. Ill CI 



In order to compute the free energy as function of 
magnetization, we consider the Hamiltonian ^ with 
an external magnetic field h, it then becomes H = 
—N{mP + hm). Applying the integral representation of 
the delta function, one finds 



Z{h) =^ f dmS{Nm~ M{{S}))e 



dm dX e/3^(™'' + '»"i)-JVAm+Ar log (2 cosh A) 



The saddle point condition imposes that A = fipmP 
/3h and therefore 



Z(h)= / dTOe-'^^^^"^'^-'^) 



(17) 



with 



g{m, /3) = (p - 1)toP - - log [2 cosh ((3pmP-'^ + I3h)] 

P 



(18) 



3. Solving the dynamics of the melting process 

We define the Glauber dynamics as follow: At each 
interval of time 1/iV, we pick up one spin, and set it plus 
with probability p+ = (1 +tanh/3/i)/2 and minus with 
probability p_ = (1 — tanh/3/i)/2, where h = pmP~^ is 
the local field on the spin, that is conveniently the same 
for every spin. To compute the evolution of magneti- 
zation and susceptibility ([7]), we define P{t,M) as the 
probability (over different realizations of the dynamics) 
that the value of the magnetization is M at time t, where 
—TV < M < N . Its averaged value per site reads 



1 

m{t) j Pit, M)Mdm , 



(23) 



and the average susceptibility is 

X{t) = 4 W M)M^dm - [ P{t, M)Mdm 

^ Jo VJo 



(24) 
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where M = Nm. The master equation for P{t, M) reads 

P{t + dt, M) = ^ + + p- {M + 2) P{t, M + 2) 
+ ^-lr!-^^ P.(M-2)Pft,M-2) 



2iV 
iV + M 



N - M 
p^{M) + ——p^{M) 



P{t,M). (25) 



27V ' 2N 

The average magnetization at time t + dt is thus 



Nm{t+dt) 



[ P{t,M)dm[M -m + ta.nh{PpmP~^)] . 
Jo 

(26) 

Considering that fluctuations in m are 0(1 /Vn), this 
leads to the differential equation (|4]). 

The evolution of the dynamical magnetic susceptibility 
is computed in a similar manner. Using (j25p . the second 
moment of M reads 

[ M^P{t + dt,M)dm^ [ M^P{t,M)dm 
Jo Jo 

+ [ 2MP{t, M) [-m + tanh {l3pmF'^)\ dm 
Jo 

- I 2P(i,Af) [mtanh(/3pmP"i) - 1] dm. (27) 



Now from the definition of x(i) after collection all the 
terms from (j26D and (l?7l) we get 



N x{t + dt)/p^Nx{t)/(3 

+ I 2{M - {M))P{t, M) [-TO + tanh (/Jprn^-^)] dTO 
Jo 

+ 2{miimh{l3pmP'^)) - 2 - (-m + tanh (/JptoP^I))^ . 



The remaining integral on the r.h.s. can be rewritten 
using a Taylor series around {m(t)). Finally we obtain 
the following differential equation for x{t) 



dt 



{cosh [/3pm(i)P-i]}2 



+ 213- I3m^{t)- p{i&n\v[Ppm{ty-^]Y . (28) 



Equations Q and (1^51) can be numerically integrated to- 
gether, and yield the results in Fig. [51 
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